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SUMMARY 


Three basic types of symmetry (and their combinations) exhibited by tire 
response are identified. A simple and efficient computational strategy is pre- 
sented for reducing both the size of the model and the cost of the analysis of 
tires in the presence of symmetry-breaking conditions (e.g., unsymmetry of the tire 
material, geometry, and/or loading). The strategy is based on approximation of the 
unsymmetric response of the tire with a linear combination symmetric and anti- 
symmetric global approximation vectors [or mode 6 ). 

The three main elements of the computational strategy are as follows: (1) use 

of three-field mixed finite-element models having independent shape functions for 
stress resultants, strain components, and generalized displacements, with the stress 
resultants and the strain components allowed to be discontinuous at interelement 
boundaries; (2) use of operator splitting (additive decomposition of some of the 
matrices and vectors in the finite-element model) to delineate the symmetric and 
antisymmetric contributions to the response; and (3) successive use of the finite- 
element method and the classic Rayleigh-Ritz technique to substantially reduce the 
number of degrees of freedom. The finite-element method is first used to generate 
a few global approximation vectors (or modes) . Then the amplitudes of these modes 
are computed with the Rayleigh-Ritz technique. 

The proposed computational strategy is applied to three quasi-symmetric prob- 
lems of tires, namely, (1) linear analysis of anisotropic tires through use of 
semianalytic finite elements, (2) nonlinear analysis of anisotropic tires through 
use of two-dimensional shell finite elements, and (3) nonlinear analysis of ortho- 
tropic tires subjected to unsymmetric loading. In the first two applications, the 
anisotropy (nonorthotropy) of the tire is the source of the symmetry breaking; in 
the third application, the quasi -symmetry is due to the unsymmetry of the loading. 
The effectiveness of the proposed computational strategy is also demonstrated with 
numerical examples, and its potential for handling practical tire problems is 
outlined. 


1 INTRODUCTION 

One of the most challenging applications of coirputational structural mechanics 
is the numerical simulation of the response of aircraft tires during taxi, take-off, 
and landing operations. In addition to the harsh environment tires are subjected to, 
their composite structure is composed of rubber and textile constituents which ex- 
hibit anisotropic, nonhomogeneous material properties. (See fig. 1.) Tires are 
subjected to inflation pressure and to a variety of unsymmetric mechanical and 
thermal loads which can result in large structural rotations and deformations as 
well as in variations in the characteristics of the tire constituents. Also, the 
laminated carcass of aircraft tires is thick enough to allow significant transverse- 
shear deformations. Moreover, the detailed stress and temperature distributions in 
tires may require the use of three-dimensional finite elements to model certain 
regions of the tire. To date only a few analytical tools have been available to 
assist the tire designer. A review of the commonly used models for predicting tire 
response is given in references 1 to 3. 

The aforementioned difficulties can make the cost of the numerical simulation 
of the tire response prohibitive. Hence, there is a need to develop modeling 
strategies and analysis methods to reduce this expense. In recent years nonlinear 
analyses of static and dynamic problems have become the focus of extensive research 


efforts. This endeavor has prompted the development of versatile and powerful finite- 
element discretization methods as well as the development of improved numerical 
methods and software systems for nonlinear static and dynamic analyses of structures 
and solids. Novel techniques which have emerged from these efforts include three- 
field mixed finite-element models (refs. 4 and 5) , semianalytic finite-element 
models for shells of revolution (refs. 6 and 7) , techniques for exploiting the 
special symmetries exhibited by anisotropic plates and shells in their finite- 
element analyses (refs. 8 and 9), single- and multiple-parameter reduction techniques 
(refs. 10 to 13)/ and operator splitting techniques (refs. 14 to 16). Current tire 
modeling studies at NASA Langley Research Center are aimed at developing an accurate 
and cost-effective computational strategy for predicting tire response. This is 
accomplished through combining the aforementioned analysis techniques with a modeling 
strategy for exploiting the symmetries exhibited by the tire response. In situations 
for which the symmetry is broken because of unsymmetry of the tire material, geometry, 
or loading (quasi-symmetric problems) , the proposed strategy allows the reduction of 
the size of the analysis model to that of the corresponding symmetric tire problem. 

The objectives of the present paper are the following: (1) to review the dif- 

ferent types of symmetry exhibited by the tire response; (2) to present a simple and 
efficient computational strategy for the analysis of tires in the presence of 
symmetry-breaking conditions; and (3) to discuss the potential of the proposed 
computational strategy and its application to practical tire problems. To sharpen 
the focus of the study, discussion is limited to tires with elliptic cross sections 
and linear material response. The tire is modeled through use of the moderate- 
rotation, geometrically nonlinear Sanders-Budiansky shell theory with the effects 
of transverse-shear deformation and laminated anisotropic material response included 
(refs. 17 and 18). The sign convention for the external loading, generalized dis- 
placements, and stress resultants is given in figure 2. 


2 SYMMETRIES EXHIBITED BY TIRE RESPONSE 

The fundamental definitions of symmetry, symmetry elements, and symmetry 
operations are reviewed in reference 19. When applied to tires, the "axiom of 
symmetry" introduced in references 19 and 20 becomes the following: given a tire 

exhibiting certain types of symmetry and a system of loads which exhibits the same 
types of symmetry as those of the tire, the response obtained is expected to exhibit 

the same types of symmetry as those inherent in the tire and the loading system. 

Herein the symmetry of the tire refers to the symmetry of (1) tire geometry, (2) 
lamination and material parameters (e.g., cord orientation and stacking of layers), 
and (3) boundary conditions at the rim. 

The symmetry of tires is described by giving the set of all operations which 

preserve the distance between pairs of points of the tire and which take the tire 

into an equivalent configuration (i.e., into a configuration which is indistinguish- 
able from the original configuration but not necessarily identical to that 
configuration). Any such operation is called a symmetry transformation. The set 
of all symmetry transformations forms the symmetry group of the tire (ref. 20). 

The three basic types of symmetry frequently exhibited by the tire response 
are (1) axial symmetry, (2) rotational symmetry, and (3) reflection symmetry. The 
characteristics of these symmetries are examined subsequently. 
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2 . 1 Axial Symmetry 


A tire is said to exhibit axial symmetry about its axis of revolution (y-axis 
in fig. 3) if it can be brought into an equivalent configuration by any rotation 
about that axis. Axial symmetry is exhibited by tires whose geometry, material 
properties, and boundary conditions are independent of the circumferential coordinate 
0 (i.e., the tire is axisymmetric) . If the loading is also axisymmetric , then the 

response is expected to be axisymmetric. (This excludes the case of unsymmetric 
buckling of a tire subjected to axisymmetric loading.) For isotropic and ortho- 
tropic tires subjected to axisymmetric loading, the two sets of generalized displace- 
ments, strain components, and stress resultants (u, w, <f> ; e , e_, K n , 2e 

s s 0 0 s3 

and N s , H e , M s , M„, S s and v, 2z^. 2*^, 2^, «nd N^, M >e> 5,) 

are generally uncoupled. In contradistinction, for anisotropic tires (also sub- 
jected to axisymmetric loading) the two sets are coupled. (See fig. 4.) In the 
presence of axial symmetry only one meridian of the tire needs to be analyzed. 

(See fig. 4.) 


2 . 2 Rotational Symmetry 

A tire exhibits rotational symmetry about an axis normal to the axis of revolu- 
tion (e.g., x-axis or z-axis in fig. 3) if a 180° rotation around this axis brings 
the tire into an equivalent configuration. Such rotations are referred to as 
dihedral rotations, in contradistinction to the axial rotation discussed in 
section 2.1. If the external loading exhibits rotational symmetry about the same 
axis as the tire, then the response will also exhibit the same symmetry. For con- 
venience, the reference (middle) surface of the tire is parametrized by the two 
coordinates 5 and 0. (See fig. 3.) The coordinate 5 represents a normalized 
arc length distance measured from the crown of the tire, and the coordinate 0 is 
the circumferential angle measured from the intersection of the negative z-axis with 
the reference surface. The rotational symmetry of the external loading and the 
response about an axis normal to the tire through the origin (5=9=0, the z-axis 
in fig. 3) are expressed by the following relations: 

External loading: 
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Generalized displacements : 



Stress resultants: 



where p , Pg , and p represent the intensities of the loading in the s- , 0 -, and 
x^-direclions to the reference surface and and are the intensities of the 

applied moments. b 

The symmetry relations for the strain components are the same as those for the 
corresponding stress resultants. These relations show that for a symmetric loading 
(defined by eqs. ( 1 )), the generalized displacements u, v, <p , and <p^ , the 
transverse-shear stress resultants Q s and Qg , and the transverse-shear strain 
components 2e s3 and 2 e 0 3 vanish at the axis of rotation (5 = 0 = 0 ). 

Rotational symmetry can be exhibited by isotropic, orthotropic, and anisotropic 
tires. It is sometimes referred to in the literature as nvCAA^ton byimojOiy because 
of its effect on the (5,0) space. The intersections of the axis of rotational 
symmetry with the tire middle surface correspond to the d&VltQJL 6 AymmoASiy in the 

( 5 , 0 ) space. 

An illustration of rotational symmetry for a two-layered tire subjected to 
uniform inflation pressure and a localized loading is given in figure 5(a), in 
which normalized contour plots for the two displacements v and w are shown. 

As can be seen from figure 5(a), only one-half of the tire needs to be analyzed. 


2 . 3 Reflection Symmetry 

A tire exhibits reflection (or mirror) symmetry with respect to a. plane (e.g.. 
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x-z plane in fig. 3) if it can be brought into an equivalent configuration by mirror 
reflection in that plane. For reflection symmetry with respect to the plane £ = o 


following: 
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The symmetry relations for the strain components are the same as those for the 
corresponding stress resultants. Reflection symmetry can occur in isotropic and 
ortho tropic tires but not in anisotropic tires. 


2 . 4 Symmetry Combinations 

Combinations of the foregoing three basic types of symmetry are possible in 
tire response. These combinations reduce the size of the analysis model and the 
associated computational cost beyond what is possible with only one symmetry 
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transformation. The symmetry combinations can be of the same type or of different 
types. The following three examples of symmetry combinations are discussed herein. 

2.4.1 Reflection symmetry combined with dihedral rotational symmetry .- This 
results in reflection symmetry with respect to two planes , as shown in figure 5(b) 
for a two-layered orthotropic tire subjected to uniform inflation pressure and 
localized loading. The size of the analysis model is only one quadrant of the tire. 

2.4.2 Periodic (or translational) symmetry .- Periodic symmetry is rotational 
symmetry around both the x- and the z-axis (see fig. 3) and is exhibited by the 
response of orthotropic and anisotropic tires, for which the tire geometry, material, 
boundary conditions, and loading have rotational symmetry around both axes. An 
example of periodic symmetry for two-layered orthotropic tires subjected to uniform 
inflation pressure and "symmetric" localized loading is given in figure 5(c). The 
size of the analysis model is only one octant of the tire. 

2.4.3 Axial symmetry combined with rotational and/or reflection symmetry .- If 
both the tire and the loading possess rotational and/or reflection symmetry in addi- 
tion to axial symmetry, then the size of the analysis model reduces to half the 
meridian. The symmetry conditions in these cases are listed in table 1. 


2.5 Linear Response of Tires Subjected to Unsymmetric Loading 

The aforementioned symmetry concepts can be used in predicting the linear re- 
sponse of tires subjected to unsymmetric loading. This is accomplished through the 
decomposition of the loading into various symmetric and antisymmetric components. 

The tire response is then obtained through summation of the responses to each of the 
symmetric and antisymmetric loading components. These responses have a pKloftl 
known symmetries (and/or antisymmetries) , and therefore only a small portion of the 
tire needs to be analyzed in each case, with a consequent large reduction in the 
computational effort. 

Note that a unique. conxeApondence ext*t* between the Aymmetntc and antl*ymmetntc 
component* o l the loading and the cosine* ponding Ae*pon*e component* . This is based 
on the principle of superposition, which is generally valid only for linear problems. 


3 SYMMETRY-BREAKING (QUASI- SYMMETRY) CONDITIONS 

Symmetry-breaking, or quasi -symmetry, conditions are \ised in the present study 
for situations in which the tire response cannot be exactly computed a* a finite Aum 
o£ A ymmetAlc and anti* ymmetsilc mode, 6. The major sources of symmetry breaking are 
unsymmetry of the tire material, geometry, boundary conditions, and loading. In 
the mechanics literature the term "asymmetry" is often used to denote unsymmetry. 
Material anisotropy is a source of reflection unsymmetry. The unsymmetry in geometry 
can be caused by the presence of unsymmetric imperfections in the tire. The contact 
of the tire with the pavement can be considered as unsymmetric boundary conditions. 
Unsymmetric loading is a source of symmetry breaking when the nonlinear tire response 
is considered. In a practical situation, a combination of the different types of 
unsymmetry may exist. 

If the unsymmetry of the tire material, geometry, or loading is small, the 
unsymmetric tire response can be obtained as a small perturbation from the corre- 
sponding symmetric response. As an example of this, for "slightly anisotropic" tires 
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the response can be obtained as a small perturbation from the orthotropic response 
(see refs. 21 and 22) , thereby reducing the size of the analysis model of the aniso- 
tropic tire to that of the corresponding ortho tropic tire. On the other hand, when 
the unsymmetry is large (e.g. , for strongly anisotropic tires), the classic pertur- 
bation technique does not work. In section 4, a computational strategy is presented 
for reducing the size of the analysis model for quasi-symmetric problems of tires to 
that of the corresponding symmetric problems. 


4 PROPOSED COMPUTATIONAL STRATEGY FOR OUASI-SYMMETRIC PROBLEMS 

In this section a computational strategy is presented for the efficient analysis 
of tires in the presence of symmetry-breaking (quasi-symmetric) conditions. The 
proposed strategy is based on approximating the unSymmetric response of the tire 
with a Linear combination o£ symmetric and antisymmetric global approximation vectors 
(or modes). The three main elements of the computational strategy are as follows: 

(1) use of three-field mixed models having independent shape functions for stress 
resultants, strain components, and generalized displacements, with the stress re- 
sultants and strain components allowed to be discontinuous at interelement boundaries; 

(2) use of operator splitting (additive deconposition of some of the matrices and 
vectors used in the finite-element model) to delineate the symmetric and antisymmetric 
contributions to the response; and (3) successive use of the finite-element method 
and the classic Rayleigh-Ritz technique to substantially reduce the number of degrees 
of freedom. The finite-element method is first used to generate <X global approxi- 
mation vectors (or modes) . Then the amplitudes of these inodes are computed with 

the Rayleigh-Ritz technique. 

Main elements 2 and 3 are essential for the success of the computational strategy, 
and main element 1 significantly enhances the efficiency and effectiveness of the 
strategy. In sections 5 to 7 the details of application of the proposed strategy 
to the following three quasi-symmetric problems are discussed: linear analysis of 

anisotropic tires (subjected to arbitrary loading) through use of semianalytic 
finite elements; nonlinear analysis of anisotropic tires through use of two- 
dimensional shell finite elements; and nonlinear analysis of orthotropic tires 
subjected to unsymmetric loading. In the first two applications the anisotropy 
(nonorthotropy) of the tire is the source of symmetry breaking; in the third appli- 
cation, the symmetry breaking is due to the unsymmetry of the loading. Numerical 
results are presented in section 8. 

The three applications presented in sections 5 to 7 are intended to demonstrate 
the use of the proposed computational strategy in analyzing practical tire problems. 
Tires with a small number of plies (e.g., two) generally exhibit a high degree of 
anisotropy. As the number of plies in the tire increases, the effect of the aniso- 
tropic (nonorthotropic) stiffness coefficients on the tire response decreases, and 
for 10 or more plies the tire can be treated as orthotropic. 

For tires with uniform circumferential properties subjected to pressure load- 
ing, the response can be accurately predicted with the linear theory. If the 
loading is nonuniform in the circumferential direction, then an efficient analysis 
technique can be obtained from a Fourier series representation of the different tire 
quantities in the circumferential direction and a one-dimensional finite-element 
discretization in the meridional direction. However, for anisotropic tires this 
analysis technique can become expensive. The application of the aforementioned 
computational strategy significantly reduces the cost of analyzing highly anisotropic 
tires, as described in section 5. 
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When the tires are subjected to a combination of pressure loading and localized 
loading (e.g., in the contact region), then the use of two-dimensional shell finite 
elements is appropriate. However, this technique can also become quite expensive 
for anisotropic tires or for unsymmetric localized loading. The use of the computa- 
tional strategy to reduce the cost of the analysis in these cases is discussed in 
sections 6 and 7. 

The computational strategy is based on a form of the moderate-rotation Sander s- 
Budiansky shell theory with the effects of transverse-shear deformation and laminated, 
anisotropic material response included. (See refs. 17 and 18 and appendix A.) A 
total Lagrangian formulation is used and the fundamental unknowns consist of the 
five generalized displacements, the eight stress resultants, and the corresponding 
eight strain components of the middle surface. (See fig. 2 for sign convention.) 


5 LINEAR ANALYSIS OF ANISOTROPIC TIRES THROUGH USE 
OF SEMI ANALYTIC FINITE ELEMENTS 


For tires with uniform circumferential properties, an approach generally used 
for predicting the linear response is based on the representation of the shell vari- 
ables and the external loadings with Fourier series in the circumferential coordinate 
0 combined with the use of finite elements in the meridional direction. (See 
ref. 23.) The following expressions are used for the external loadings: 


P s (s,0) 
P 0 (s,0) 
p(s,0) 


m (s,0) 
s 

m (s,0) 


oo 

- 


cos n0 


oo 

♦ 2 > 


Pe,n 


sin n0 


The generalized displacements X, stress resultants H, and strain components E 
are approximated with expressions of the form 


oo w 

5,0) = / . W i X i cos n0 + / ^ N^X^ 

“ n *—r n 


sin n0 


oo 00 

[ (s r 6) = / . MV" cos n0 + 

n n 


sin n0 


oo oo 

E (s , 0) = E cos n0 + z NV- 

r\ ^ n n 


sin n0 


where W 1 are the shape functions used in approximating the generalized displace- 
ments in the meridional direction; Sf 1 are the shape functions used in approximating 

the stress resultants and the strain components? X 1 and X 1 refer to the gener- 

n n si 

alized nodal displacement coefficients associated with the Fourier harmonic n; H 
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and H n refer to the stress-resultant parameters associated with the harmonic n; 

and and refer to the strain parameters associated with the harmonic n. 

n n 

The barred tire quantities are associated with the sine series. Based on the pres- 
ent study, the degree of polynomial shape functions is one lower than that 

of N 1 . Moreover, the continuity of the stress resultants and the strain components 
is not imposed at the interelement boundaries and, therefore, both the strain param- 
eters and the stress-resultant parameters can be eliminated at the element level. 

In equations (7) to (10) , the range of the superscript i is 1 to m, the 
number of displacement nodes in the element; the range of superscript Z is 1 to 
6, the number of parameters used in approximating each of the stress resultants or 
the strain components. A repeated superscript denotes summation over its entire 
range. 


5.1 Governing Finite-Element Equations 

The governing finite-element equations are obtained through application of the 
three-field Hu-Washizu mixed variational principle. For tires with uniform circum- 
ferential properties, the finite-element equations uncouple in harmonics, that is, 
each load harmonic excites only the response in the same harmonic. For the harmonic 
n, the governing equations for individual elements can be cast in the following form: 



where {x} , {e} , and {h} are respectively the vectors of nodal displacement 

n n n 

coefficients, the strain parameters, and the stress-resultant parameters associated 

with n; [k] , [r] , and [s] are "generalized" stiffness matrices; {Q> is 

n n n n 

the vector of normal external loadings; subscript n refers to the nth Fourier 

harmonic; and superscript t denotes transposition. 


If the strain parameters are eliminated from equations (11) on the element 
level, one obtains the governing equations of the two-field Hellinger-Reissner 
mixed models. The governing equations for the individual elements of these models 
take the form 



where [f] is the linear flexibility matrix of the element given by 

[f] - WhKrtR] 

The equivalence between the three-field Hu-Washizu mixed models, the two-field 
Hellinger-Reissner mixed models, and the reduced-selective integration displace 
ment models is discussed in references 4 and 24. 
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5 . 2 Operator Splitting 


To simplify the analysis of anisotropic tires, the vectors {E} n , {H} n , 

{x} n , and (Q} n in equations (11) are each partitioned into two subvectors : 

{£} , (5) ; (H) , {H} ; (X) / {X} ; and {{}} , {Q.} . The first and second sub- 

nnnnnn nn 

vectors of each pair are associated with the symmetric and antisymmetric fundamental 

unknowns (or external loadings) , respectively. The two sets of symmetric and anti- 
symmetric shell variables are listed in table 2. In table 2, the shell variables 
without a bar are the coefficients of the cosine series, and the barred tire vari- 
ables are the coefficients of the sine series. (See eqs. (8), (9), and (10).) 

The matrix [k]^ equations (11) is now partitioned into two block-diagonal 

submatrices [K]^ and associated with {£ and {E} r , respectively, and 

a coupling submatrix [k]^. The coupling submatrix contains all the nonorthotropic 

material stiffnesses of the tire and, therefore, vanishes for orthotropic (and 
isotropic) tires. 

The matrices [ R l n and [s] n are partitioned into the submatrices 

[R] n , [RL and [#] , [5] associated with {£} , {H} , {X} and {€} , {H} , 

nn nn _nnn nn 

(X} n . No coupling submatrices other than [lt] n exist between the two sets. 

The governing finite-element equations for each individual element are parti- 
tioned into two sets of coupled equations and are embedded into a single-parameter 
family of equations of the form 




(13) 


(14) 


In equations (13) and (14) , X is a tracing parameter which identifies all the non- 
orthotropic (anisotropic) contributions contained in the matrix [k] r . For X = 0, 
equations (13) and (14) are uncoupled; that is, the response to symmetric loading 
is uncoupled from the response to antisymmetric loading. 

The explicit forms of the arrays appearing in equations (13) and (14) are given 
in appendix B. For n > 0, [K] n = [K] n , [R] n = [R] n , [5] n = [S] n , and 

ra n -R£- 


5.3 Application of Reduction Method 

A reduction method is now applied to uncouple governing equations (13) and (14) 
and to reduce the total number of degrees of freedom of the initial discretization. 
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The method is based on generating a few global approximation vectors or modes at 
A = 0 with the finite model- The vectors of nodal displacement coefficients , 
stress-resultant parameters, and strain parameters are then expressed as linear 
combinations of these global approximation vectors in the following transformation: 


W (15) 

n 

n 

and 


W (16) 

n 
n 

where the matrices [r] are transformation matrices and {tD } is a vector of 

n r n 

undetermined coefficients (amplitudes of global approximation vectors) which are 
functions of A. 




An effective choice for the global approximation vectors (the columns of the 
matrices [r]^) consists of the solution corresponding to A = 0 and its various- 

order derivatives with respect to A (evaluated at X = 0) , that is. 







(E, „ ( I } n 


(H) {|£} 

n dX n 


{ X} 

n 3A n 



(17) 


(18) 

(19) 


with similar expressions for [iV] , [r„] > and [r y ] in terms of {E}, {H}, 

£ n ft n An 

and {X}. The total number of global approximation vectors is K. 


The equations used in generating the global approximation vectors are obtained 
by successive differentiation of the finite-element equations of the discretized 
shell (eqs. (13) and (14)) and are listed in appendix C. The global approximation 
vectors are generated at A = 0, thereby uncoupling equations (13) and (14) . 

For symmetric loadings the vectors {E}, {H}, {X}, and all their even- 

order derivatives with respect to A are zero- Also, the odd-order derivatives 
of {£} , {H}, and {X} with respect to A vanish. For antisymmetric loadings, 

the situation is reversed with respect to the derivatives of {E} , {H} , and {X} 

and of {£}, {H}, and {X}. 
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The classic Rayleigh-Ritz technique is used to replace the governing finite- 
element equations of the shell with the following reduced system of equations in 
the unknown parameters (components of the vector { 1 ^}^) : 

[[k]„ ♦ >M„] (*}„ = (20) 

where 

[k] - E ([r E ]‘W[r t ] + [r r ] t [R][r„] 

Elements t n 

+ [r^Mtrg] + [r E ] t [R][r (( ] + Cr J( D t [a] t Cr E ] + Cr H ] t [»][r J( 3 
+ Cr H ] t [R] t [r E ] + Cr (f 3 t [s3[r x D + [r Jf ] t [$] t [r I( ] + Cr x 3 t ts] t tr H 3) (2D 

[t] - E (Cr^raCrj] + [r E 3 t K3 t Cr _]) (22) 

Elements 

and 

= Y (Cr x ] t {<9} + [r x ] fc { 0 ) (23) 

Elements 

In equations (21) to (2 3) the subscript n has been dropped for convenience. 

Note that the chosen global approximation vectors (eqs. (17) to (19)) are the 
same as those used in the Taylor series expansion of the response in terms of X. 
These vectors provide a direct measure of the sensitivity of the different response 
quantities to nonorthotropic (anisotropic) material stiffness coefficients of the 
tire. A comparison between the predictions of the foregoing technique and the 
Taylor series e^qpansion is given in section 8. 


5.4 Case of Axi symmetric Loading 

If the Fourier harmonic n is set equal to zero in equations (7) to (10) , the 
governing equations for the axisymmetric loading case are obtained. All the barred 
loading components and shell variables vanish, and the subvectors {£}, {E}, {H}, 

{ H} r {X}, {X}, {<§}, and {Q} are now associated with the fundamental unknowns 

(and external loadings) , as shown in table 3. 

The corresponding matrices in equations (13) and (14) have different dimensions. 
For X = 0, the uncoupling of equations (13) and (14) signifies the uncoupling of 
the response of orthotropic (and isotropic) shells to normal and meridional loading 
from the response to circumferential loading. (See sections 2.1 and 2.4.) 


6 NONLINEAR ANALYSIS OF ANISOTROPIC TIRES THROUGH USE 
OF TWO-DIMENSIONAL SHELL FINITE ELEMENTS 

As a second application of the aforementioned computational strategy, consider 
the geometrically nonlinear analysis of anisotropic tires by use of three-field 
mixed two-dimensional shell finite elements. The material stiffness matrix of the 
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shell is decomposed into the sum of an orthotropic and a nonortho tropic (anisotropic) 
part, and the governing nonlinear finite-element equations for individual elements 
are embedded into a two-parameter family of equations of the form 


( 


K 

R 1 

0 


R 

0 


0 

s 

0 


K 

0 

0 



(24) 


where {x}, {e}, and {h} are the vectors of nodal displacements, strain param- 
eters, and stress-resultant parameters; (G(X)} and {flt(H,X)} are vectors of 
nonlinear contributions; {q} is the vector of normalized external loading 
coefficients; q is a loading parameter; and A is a tracing parameter which 
identifies all the nonorthotropic (anisotropic) ^contributions to the element 
equations. The explicit forms of the arrays [K], [K]> [rL [s], (G(x) }, 

(flt(H,X) }, and {q} are given in appendix D. 


With A and q chosen as the control parameters, a two-parameter reduction 
method is applied to reduce both the size of the analysis model and the total 
number of degrees of freedom used in the initial discretization. This reduction 
is accomplished by use of the following transformation; 



(25) 


where is a vector of amplitudes of global approximation vectors which are 

functions of both A and q, and the [r] terms are transformation matrices 
whose columns are the global approximation vectors consisting of the solution at 
A = q = 0 and its various-order derivatives with respect to A and q, that is, 



= [{e} 

{ M, 

{M} 

'ax' 

*3q ax> 

...] 

(26) 


= [{H} 

{ M } 

{Si} 

'ax' 

<4 


(27) 


= [{X} 

{Si, 

(M, 

ax 

<arfx> 

...] 

(28) 


The equations used in generating the global approximation vectors are obtained 
by successive differentiation with respect to q and A of the finite-element 
equations of the discretized tire (eqs. (24)) and are listed in appendix E. For 
convenience, the initial set of global approximation vectors is generated at 
A = q = 0. This reduces the computational effort in evaluating the global approxi- 
mation vectors because of the following: for q = 0, {e}, (h } , and {x} vanish, 

as do all the nonlinear terms on the left-hand sides of equations (El) ; and for 
A = 0, each of the global approximation vectors exhibits known types of symmetries 
(or antisymmetries) similar to those exhibited by the response of orthotropic tires, 
and therefore the size of the model used in generating these vectors is the same as 
that used for orthotropic tires. 
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The classic Rayleigh-Ritz technique is used to replace equations (24) with 
the following system of nonlinear equations in the unknown parameters { 1 ( 1 }: 


Jjk] + x[fe]J{>} + (m(ij>)} - q{q} = o 

where 

[it] - E (tr E ] t [K][r E ] + [r ]‘[R][r ] 

Elements 

* + [r H ] t [s][r x ] + [r x ] t [s] t [r H ]) 

M - E 

Elements 


( 29 ) 


(30) 

(31) 


{mw> = ( [r H ] t {©(^) } + [r x ] t {l!l(4>)}) 

Elements 


(32) 


and 

2 t r x ^{Q> (33) 

Elements 

The nonlinear vectors {G(^)} and {flt(^) } are obtained from the vectors (G(X)} 
and {ffl(H,X) } through replacement of {h} and {x} with their expressions in 
terms of {x^ } (eqs. (25)). 


7 NONLINEAR ANALYSIS OF ORTHOTROPIC TIRES 
SUBJECTED TO UNSYMMETRIC LOADING 

As a final application of the computational strategy, consider the geometri- 
cally nonlinear analysis of orthotropic tires subjected to the localized loading 
shown in figure 5(b). The loading is unsymmetric with respect to reflection in 
the x-y plane, but it is symmetric with respect to reflections in the x-z and y-z 
planes. (See fig. 3.) The tire is modeled by use of three-field mixed two- 
dimensional shell elements. The loading vector is decomposed into symmetric and 
antisymmetric loadings (with respect to reflections in the x-y plane) , and the 
governing nonlinear finite-element equations for individual elements are embedded 
in a two-parameter family of equations of the form 



where q^ and q 2 are normalized loading parameters and {Q^} and are the 

vectors of normalized symmetric and antisymmetric loadings. 

With and q^ chosen as the control parameters, a two-parameter reduction 

method is applied to reduce both the size of the analysis model and the total 
number of degrees of freedom used in the initial discretization. This reduction 
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is accomplished with the procedure outlined in section 6. Note that the evaluation 
of the global approximation vectors at q 1 = = 0 allows the reduction of the 

size of the analysis model to just one octant of the tire (instead of the one 
quadrant required for the original unsymmetric loading) . Each of the global ap- 
proximation vectors exhibits known symmetry (or antisymmetry) along the centerlines. 

The following observations can be made with regard to the proposed computational 
procedure : 

1. Although the individual global approximation vectors exhibit symmetry (or 
antisymmetry) along their centerlines, their linear combination, which approximates 
the response to the general unsymmetric loading, is generally unsymmetric. 

2. The reduced equations can be formed by use of only one octant of the tire. 
However, the solution of these equations approximates the unsymmetric response over 
one quadrant of the tire and consequently over the entire tire. 

3. In case the given loading is antisymmetric (i.e., q^ = 0) , the procedure 
can be simplified, since only one parameter q is needed. Obviously, for symmetric 
loadings the response is symmetric. 

4. Updated sets of global approximation vectors may be formed with the compu- 
tational strategy, but it is essential that the tire configuration used in generating 
these vectors be symmetric (e.g. , q^ = 0) . 

8 NUMERICAL EXAMPLES 

To evaluate the effectiveness of the proposed computational strategy, a number 
of quasi-symmetric tire problems have been analyzed. For each problem, the solutions 
obtained with the proposed strategy are compared with the direct analysis of the 
tire. The results of three typical problems are discussed herein. The three 
problems are linear response of an anisotropic tire subjected to normal pressure 
loading, linear response of an anisotropic tire subjected to uniform inflation pres- 
sure and localized loading, and geometrically nonlinear response of an orthotropic 
tire subjected to uniform inflation pressure and unsymmetric localized loading. The 
three problems correspond to the computational procedures presented in sections 5 
to 7. 


In the first problem, the tire is discretized with semianalytic finite elements 
(Fourier series expansions in the circumferential direction and finite elements in 
the meridional direction). In the last two problems, two-dimensional shell finite 
elements are used for the discretization. The tires are assumed to be rigidly 
clamped at the rim (at £ = ±1), to have elliptic cross sections, and to be made 
of cord-rubber composite material. The geometric and material properties are given 
in figure 6. Small research-oriented programs were used to obtain all the numeri- 
cal results presented herein. 


8.1 Linear Response of an Anisotropic Tire Subjected 
to Normal Pressure Loading 

The first problem considered is that of the linear response of an anisotropic 
tire subjected to normal pressure loading which is uniform in the meridional 
direction. Two loading cases are considered. The first is an axisymmetric loading 
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(inflation pressure) with n = 0, and the second is a nonaxi symmetric loading with 
n = 2. Since both the tire and the loading have axial symmetry as well as dihedral 
rotational symmetry, only half the meridian is analyzed. The symmetry conditions 
at 5 = 0 are u = $ =0 and v = w = 4 >q = 0. Eight one-dimensional finite 

elements are used to obtain converged solution. Quadratic interpolation functions 
are used for approximating each of the stress resultants and strain components, and 
cubic interpolation functions are used for approximating each of the generalized 
displacements. The nonzero degrees of freedom consist of 118 generalized displace- 
ment coefficients, 192 stress parameters, and 192 strain parameters. Typical 
results are presented in figures 7 and 8 for n = 0 and in figure 9 for n = 2. 

The proposed strategy is applied to this problem and the global approximation 
vectors are generated at X = 0 and used to generate the solutions at X = 1. 

The distribution of the first three nonzero global approximation vectors along 
the meridian for n = 0 is shown in figure 7. Note that the odd-order derivatives 
of {£}, {X}, and {H} and the even-order derivatives of {E}, {X}, and {H} 

with respect to X are zero. 

The even-order derivatives of u , w_, N , and M are shown in 

u o y,o y,o 

figure 7(a), and the odd-order derivatives of v^, <f>„ N A and M A ^ are 

0 U , 0 so , 0 So,0 

shown in figure 7 (b) . The magnitudes of the first-order and higher order de- 
rivatives of the different response quantities are an indication of the sensitivity 
of these response quantities to the anisotropic material coefficients. 

The tire responses corresponding to X = 1 and X = 0 are shown in figure 8 
for n = 0 and in figure 9 for n = 2, along with the predictions of the proposed 
strategy based on the use of four and six (or six and eight) global approximation 
vectors. The symmetric response quantities (with respect to 0=0) are shown in 
figures 8(a), 9(a), and 9(c), and the antisymmetric response quantities are shown 
in figures 8(b), 9(b), and 9(d). Note that for n = 0, all the antisymmetric re- 
sponse quantities corresponding to X = 0 vanish. 

As shown in figures 8 and 9, for both n = 0 and n = 2 the response of 
the two-layered tire is sensitive to the nonorthotropic stiffness coefficients 
identified by X. As the number of layers increases, the response becomes less 
sensitive to X. 

An indication of the accuracy of the solutions obtained with the proposed 
strategy is given in figure 8 for n = 0 and in figure 9 for n = 2. The gen- 
eralized displacements and stress resultants predicted with the proposed strategy 
through use of six vectors are highly accurate except for N^ for which 

eight vectors are needed. (See fig. 9(b).) The solutions obtained with eight 
vectors are almost indistinguishable from the direct finite-element solution. 

The numerical values of the amplitudes of the global approximation vectors 
in the proposed strategy are given in table 4. The corresponding amplitudes of 
the same vectors in the Taylor series expansion are also given in table 4. 

8.2 Linear Response of an Anisotropic Tire Subjected to Uniform 
Inflation Pressure and Localized Loading 

The second problem considered is that of the linear response of the aniso- 
tropic tire shown in figure 6 when subjected to uniform inflation pressure and 
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localized loading (simulating the contact pressure) . Typical numerical results 
are shown in figures 5(a), 10, and 11. 

In figure 5 (a) normalized contour plots are shown for the two displacement 
components v and w produced by the given loading. In figure 10, normalized 
contour plots are shown for the generalized displacements and their first four 
derivatives with respect to X evaluated at X = 0. Figure 11 shows the accuracy 
and convergence of the normal displacement at the center w and of the total 
strain energy U obtained with the proposed computational strategy. 

As expected, the normal displacement w, the stress resultants N , N„ , 

S o 

N . , M , M , and M _, and the strain components e , e„, 2e . , k , and 

SU S U SO S U S0 S y 

2 k s q exhibit rotational (inversion) symmetry. The in-plane displacements u and 

v, the rotation components <\>^ and (j)^, the transverse-shear stress resultants 

and Q , and the transverse-shear strains 2z _ and 2e_ exhibit rotational 
U S3 u 3 

(inversion) antisymmetry. These symmetry properties are used in conjunction with 

the procedure outlined in references 7 and 8 to reduce the analysis model to half 

the tire. (See fig. 5(a).) A uniform grid of two-dimensional elements is used. 

Biquadratic shape functions are used to approximate each of the generalized 

displacements, and bilinear shape functions are used to approximate each of the 

stress resultants and the strain components (a total of 717 nonzero displacement 

degrees of freedom, 1280 stress-resultant parameters, and 1280 strain parameters). 

The Gauss-Legendre quadrature formula with four quadrature points in the element 

domain is used for the evaluation of the elemental arrays. The highly anisotropic 

response of the tire is given in figure 5(a), in which normalized contour plots are 

shown for the v and w displacements. 

The proposed strategy is applied to the analysis of the tire. Global approxi- 
mation vectors are generated at A = 0 and are used to generate the solution 
corresponding to A = 1. The global approximation vectors exhibit reflection 
symmetry (or antisymmetry) with respect to £ = 0 and 0 = tt. These symmetry 
relations are shown in figure 10. They clearly demonstrate the fact that the global 
approximation vectors can each be obtained by analyzing only one quadrant of the 
tire. 


An indication of the accuracy and convergence of the solutions obtained with 
the proposed strategy and with the Taylor series expansion is given in figure 11. 

The standard of comparison is taken to be the direct linear finite-element solution 
of the anisotropic tire with the technique described in references 8 and 9. As 
can be seen from figure 11, the solutions obtained with the proposed strategy are 
considerably more accurate than those obtained with the corresponding Taylor series 
expansions (based on the same global approximation vectors). As an example of this, 
the errors in the total strain energy obtained with the Taylor series with four and 
seven global approximation vectors are 16.3 percent and 5.25 percent. The corre- 
sponding errors in the strain energy obtained with the proposed strategy are only 
0.15 percent and 0.001 percent. 

8.3 Geometrically Nonlinear Response of an Orthotropic Tire Subjected 
to Uniform Inflation Pressure and Unsymmetric Localized Loading 


The last problem considered is that of the geometrically nonlinear response 
of a tire subjected to uniform inflation pressure and localized loading which is 
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unsymmetric with respect to reflection in the x-y plane. (See fig. 3.) The maxi- 
mum intensity of the localized loading is 8.662 x lo 1 * * * 5 Pa (125 psi) . The material 
properties are given in figure 6, and to simplify the analysis, the anisotropic 
bending-extensional coupling coefficients of the tire are neglected (i.e., the 
tire is treated as orthotropic) . 

The unsymmetric localized loading is decomposed into symmetric and antisymmet- 
ric components, as shown in figure 12. The path parameters are selected to be the 
two loading parameters q^ and q 2 associated with the symmetric and antisymmetric 
loading components. Only one quadrant of the tire is modeled with a uniform grid 
of two-dimensional elements (a total of 371 nonzero displacement degrees of freedom, 
640 stress-resultant parameters, and 640 strain parameters). Ten global approxi- 
mation vectors are generated at zero values of q 1 and q 2 (but at a uniform 
inflation pressure p of 6.895 x 10 5 Pa (100 psi)). The global approximation 
vectors include the axisymmetric nonlinear response ({e}, {h}, and { X } ) evaluated 

at p Q = 6.895 x io 5 Pa and q = q 2 = 0 and all its derivatives up to the fifth 

derivatives with respect to q^ and q 2 * These global approximation vectors and 
the reduced equations are generated through use of only one quadrant of the tire 
and the appropriate reflection symmetry and antisymmetry conditions. Normalized 
contour plots for the generalized displacements are given in figure 12. The re- 
flection symmetry (or antisymmetry) for the solution at q^ = q 2 = 0 and its first, 
second, and third derivatives with respect to 5 = 0 and 6 = tt/2 is presented 
in figure 12. Since the major cost of the analysis is generating the global 
approximation vectors and the coefficients of the reduced equations, the proposed 
computational strategy enhances the efficiency of the reduction technique. 

An indication of the accuracy and the convergence of the solutions obtained 
with the proposed strategy is given in table 5. The standard of comparison is 
taken to be the direct nonlinear finite-element solution of the orthotropic tire. 

As shown in table 5, the solutions obtained with the proposed strategy are highly 
accurate. The errors in the normal displacement w c and in the total strain 
energy U obtained by use of six global approximation vectors (solution at 
q^ = q 2 = 0 and all its first and second derivatives with respect to q^ and q 2 ) 

are only 0.83 percent and 0.06 percent. These errors reduce further when more 
approximation vectors are used. (See table 5.) 


9 POTENTIAL OF PROPOSED COMPUTATIONAL STRATEGY 

The proposed strategy generates the solution of a complex tire problem (with 
anisotropic materials, unsymmetric imperfections, and unsymmetric loading) using 
small or Za/ige, peAZuAbatZo n6 from the response of a simpler system (e.g., an 
orthotropic tire with no imperfections) . The strategy has high potential for 
handling practical tire problems. In particular, the following observations and 
extensions of the strategy can be made. 

1. The tire response, in both the proposed strategy and the Taylor series 
expansion of the response in terms of A, is expressed as a linear combination of 

the same global approximation vectors. However, the coefficients of the linear 

combination in the Taylor series expansion are fixed and are equal to 1, A, 

2 r-1 

A /2!, ..., A / (r - 1) ! . By contrast, the coefficients in the proposed strategy 
are left as free parameters and are determined with the Rayleigh-Ritz technique. 
Numerical experiments have shown that the use of free parameters leads to accurate 
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solutions not only within the radius of convergence of the Taylor series but also 
well beyond it. 

2. The proposed strategy can be used to analyze tires with unsymmetric 
inperfections. The unsymmetric inperfections are decomposed into symmetric and 
antisymmetric components. The terms in the governing finite-element equations (24) 
associated with the antisymmetric imperfections are grouped together and multiplied 
by the tracing parameter A. The global approximation vectors are generated at 

A = 0. The size of the analysis model used in generating these vectors is the 
same as that for a tire with symmetric imperfections. 

3. The proposed strategy can be applied to contact problems of tires. The 
contact conditions can be treated as displacement-dependent unsymmetric boundary 
conditions. The unsymmetric response of the tire is approximated by a linear com- 
bination of symmetric and antisymmetric global approximation vectors (or modes) . 

The governing nonlinear finite-element equations are restructured to delineate the 
contributions to the symmetric and antisymmetric responses. A tracing parameter X 
is used to identify all the contributions to the antisymmetric response. The size 
of the analysis model used in generating the global approximation vectors is identi- 
cal to that of the corresponding structure with symmetric contact conditions. The 
amplitudes of the global approximation vectors are obtained through solution of a 
reduced set of nonlinear equations corresponding to X = 1 . 

4. The proposed strategy can be extended to tire problems with more than one 
source for the symmetry breaking through attachment of a tracing parameter with the 
terms in the governing finite-element equations associated with each source. How- 
ever, as the number of independent tracing parameters increases, both the number of 
global approximation vectors that need to be generated and the size of the reduced 
system of equations increase at a faster rate. Therefore, for practical applications 
to tire problems, an attempt should be made to minimize the number of independent 
tracing parameters by treating some of the tracing parameters as functions of others. 

5. In the application of the proposed strategy to tires subjected to unsym- 
metric loading, it is essential that the tire configuration used in generating the 
global approximation vectors be symmetric. This may limit the applicability of the 
strategy to cases in which the magnitudes of the unsymmetric loading components 
are not very large. 

6. The proposed strategy can be applied to other mechanics problems. In 
particular, it can be used to generate the response of a complex structural system 
as a large perturbation from the response of a simpler system. It is also possible 
to use a hierarchy of simpler systems to generate the response of the complex 
structural system by choosing a number of tracing (or control) parameters and suc- 
cessively applying a single-parameter reduction method with each of the parameters. 

10 CONCLUSIONS 

Three basic types of symmetry (and their combinations) exhibited by tire 
response are identified. An efficient computational strategy is presented for 
reducing both the size of the model and the cost of the analysis of tires in the 
presence of symmetry-breaking conditions (e.g. , unsymmetric material, geometry, or 
loading) . The strategy is based on approximating the unsymmetric response of the 
tire by a linear combination of symmetric and antisymmetric global approximation 
vectors (or modes) . 
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The three main elements of the computational strategy are the following: (1) 

use of three-field mixed models having independent shape functions for stress re- 
sultants, strain components, and generalized displacements, with the stress 
resultants and the strain components allowed to be discontinuous at interelement 
boundaries; (2) use of operator splitting (additive decomposition of some of the 
matrices and vectors used in the finite-element model) to delineate the symmetric 
and antisymmetric responses; and (3) successive use of the finite-element model 
and the classic Rayleigh-Ritz technique to substantially reduce the number of de- 
grees of freedom. The finite-element method is first used to generate a few global 
approximation vectors (or modes) . Then the amplitudes of these modes are computed 
through use of the Rayleigh-Ritz technique. 

The proposed computational strategy was applied to the following three quasi- 
symmetric problems of tires: 

1. Linear analysis of anisotropic tires (subjected to arbitrary loading) 
through use of semianalytic finite elements 

2. Nonlinear analysis of anisotropic tires through use of two-dimensional 
shell finite elements 

3. Nonlinear analysis of orthotropic tires subjected to unsymmetric loading 

In the first two applications, the anisotropy (nonorthotropy) of the tire is the 
source of the symmetry breaking; in the third application, the quasi-symmetry is 
due to the unsymmetry of the loading. 

The effectiveness of the proposed computational strategy is demonstrated with 
numerical examples, and its potential for handling practical tire problems is 
discussed. The results of the study suggest the following conclusions relative 
to the importance of the three main elements of the proposed computational strategy 
and to the potential of the proposed computational strategy: 

1. The use of the three-field mixed models enhances the effectiveness of the 
proposed computational strategy and offers the following advantages over the 
equivalent displacement models: 

A. To reduce the model size of anisotropic tires, the decomposition of the 
material stiffness matrix of the tire affects only the linear arrays [K] 
and [K] of the tire. By contrast, if a displacement model is used, both 
linear and nonlinear arrays are affected by the decomposition. 

B. The generation of the global approximation vectors is easier and involves 
fewer arithmetic operations than in the displacement formulation. This is 
because the nonlinear terms of the mixed models are bilinear (or quadratic) 
in the nodal displacements and in the stress-resultant parameters and are 
not cubic as in the displacement formulation. 

C. Experience with reduction methods based on mixed models and displacement 
formulations has demonstrated that, for a given number of global approximation 
vectors, the accuracy of the solutions obtained with the reduction method and 
mixed models is higher than that obtained with the corresponding reduction 
method and displacement models. This is particularly true for stress re- 
sultants and strain components. 
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2. The use of path derivatives as global approximation vectors leads to accu- 
rate solutions with a small number of vectors. Therefore , the time required to solve 
the reduced equations is relatively small, and the total analysis time to a first 
approximation equals the time required to evaluate the global approximation vectors 
and to generate the reduced equations. The operator splitting technique allows the 
reduction of the size of the analysis model of the tire and, therefore, the time 
needed to generate the global approximation vectors is reduced. 

3. The global approximation vectors provide a direct measure of the sensitivity 
of the different response quantities to the tracing parameters used and their associ- 
ated effects (e.g., anisotropic (nonorthotropic) material coefficients of the tire, 
geometric imperfections, or unsymmetric contact conditions). The sensitivity of the 
global response can also be assessed with these vectors. 

4. The reduction method used in the proposed computational strategy exploits 
the best elements of the finite-element method and of the Rayleigh-Ritz technique, 
as follows: 

A. The finite-element method is used as a general approach for generating 
global approximation vectors. The size of the analysis model used in gener- 
ating these vectors is the same as that of the corresponding orthotropic panel. 

B. The Rayleigh-Ritz technique is used as an efficient procedure for minimizing 
and distributing the error throughout the structure. 

5. The reduction method extends the range of applicability of the Taylor 
series expansion by relaxing the requirement of using small loading and tracing 
parameters in the expansion. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
January 7, 1987 
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APPENDIX A - FUNDAMENTAL EQUATIONS OF SHELL THEORY 
USED IN THE PRESENT STUDY 

The fundamental equations of the Sanders-Budiansky type shell of revolution 
used in the present study are summarized herein. The effects of laminated, aniso- 
tropic material response and transverse-shear deformation are included. 


Strain Displacement Relationships 

The relationships between the strains and generalized displacements of the 
middle surface are given by 
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where and are the extensional strains in the meridional and circumferential 

directions, 2€ s q the in-plane shearing strain, k and are the bending 

strains in the meridional and circumferential directions, 2 k s q is the twisting 
strain, 2 £ g 3 and 2 £q 3 are the transverse-shear strains, = 3/3s, 3^ = 3/30, 

and (f> is the rotation around the normal to the shell, which is given by 


1 

2 


- | V + (a s + ~) 


(A9) 


22 


The nonlinear terms which account for moderate rotations are underlined with dashes 
in equations (Al) to (A3) . 

Constitutive Relations 

The shell is assumed to be made of a laminated, anisotropic, linearly elastic 
material. Every point of the shell is assumed to possess a single plane of elastic 
symmetry parallel to the middle surface. The relationships between the stress 
resultants and the strain measures of the shell are given by 




(A10) 


where c, f, and d are shell stiffness coefficients. The nonorthotropic 
(anisotropic) terms are circled and the dots refer to zero terms. 
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APPENDIX B - FORMULAS FOR THE COEFFICIENTS IN THE FINITE -ELEMENT 
EQUATIONS FOR SEMIANALYTIC INDIVIDUAL ELEMENTS 

The explicit forms of the elemental arrays [K] n , [K] n , [k]^, [R] n , 

[5] # [5] » {<5).,/ and {0 , associated with the harmonic n are given in this 
n n n n 

appendix. For convenience, each of these arrays is partitioned into blocks corre- 
sponding to contributions from individual nodes, stress resultants, or strain 
approximation functions. The expressions of the typical partitions (or blocks) are 
given subsequently. Note that the order of the strain parameters in these partitions 
is e g , £-0 r ^s 9 f , 2 s^^, an d 2 e 0 ^j the order of the stress- 

resultant parameters is N g , N 0 , N s0 , M , M Q , M s0 , Q , and Q 0 ; and the order 

of the generalized displacement coefficients is u, v, w, <J> , and <j> 0 . 

A typical ^Lj partition of the elemental array [K]^ is given by 
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sin n 0 


r d0 ds (Bl) 


The dots in equations (Bl) and succeeding equations refer to zero terms. The <Lj 

partition of the elemental array [K] is given by equations (Bl) after inter- 

2 2 n 
changing cos n8 and sin n0. 
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is given by 


The Ij partition of the elemental array [k]^ 
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The partition of the elemental array [R] can be written as the sum 

of two diagonal matrices as follows: 
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The partition of the array [R] is given by equations (B3) after inter- 

J o 2 n 

changing cos z n0 and sin n6 . 

The Xj partition of the elemental array [$] n is given by 


(B2) 


(B3) 
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The l j partition of the elemental array [S] 
interchanging cos 2 n0 and sin 2 n0 . n 


is given by equations (B4) after 


The ith partition of the array {Q) n is given by: 

r j 2 ' 

p cos n0 
s f n 

p;? sin 2 n9 


n 2ir 

N i W j 

u u 


■ 0 ,n 

P cos 2 n0 \ r d0 

■ r'l / 


trP cos 2 n0 
s ,n 

m Q ,n sin 2 n0 ) 


ds 


(B5) 


The ith partition of the array {0} is given by equation (B5) after inter- 

2 2 n 
changing cos n0 and sin n0 . 

In equations (Bl) to (B5) , £ is the length of the element; and are 

the shape functions for both the strain components and the stress resultants; \P 
are the shape functions for the generalized displacements; c, f, and d are shell 
stiffnesses? r is the normal distance from the tire axis to the reference sur- 
face (see fig. 1) ; and are principal radii of curvature in the meridional 

and circumferential directions ; s and 0 are the meridional and circumferential 

coordinates of the shell; and 3 = 3/3s. 

s 
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APPENDIX C - EVALUATION OF GLOBAL APPROXIMATION VECTORS 
BASED ON SEMI ANALYTIC FINITE ELEMENTS 


The global approximation vectors are chosen to be the solutions of equations 

(7) and (8), {£} » {H} , {X} , {E} , {H} , {X} , and their various-order 
n n n n n n 

path derivatives (derivatives of these solutions with respect to X) evaluated at 
X = 0. For X = 0, equations (7) and (8) are uncoupled. The path derivatives 
are obtained by successive differentiation of the governing finite-element equa- 
tions (7) and (8) with respect to X and solution of the resulting two sets of 
uncoupled equations corresponding to X = 0. For the individual finite elements, 
the recursion formulas for the path derivatives can be written in the following 
compact form: 



and 


K R 
R t 0 
0 S t 



(C2) 


where m > 1 and the vectors {P} and {P} on the right-hand sides of equa- 
tions (Cl) and (C2) are given by 1 


a m-l 

{P> = -m [k] - — - {E} 
n J n ^m-1 n 


and 


m n - 


■m [J ] 1 3 


m-1 


n 8X m " f {£>n 


(C3) 


(C4) 


Note that equations (Cl) and (C2) are uncoupled and that the coefficient matrices 
on the left-hand sides of these equations, which must be factored for each Fourier 
harmonic n, are the same for all the path derivatives. 


28 


APPENDIX D - FORMULAS FOR THE COEFFICIENTS IN THE GOVERNING 
NONLINEAR EQUATIONS FOR INDIVIDUAL TWO-DIMENSIONAL 
SHELL FINITE ELEMENTS 


The explicit forms of the arrays [K], [r], [s], [K], (G(x) }, {flt(H,x)>, 

and {Q} are given in this appendix. For convenience, each of these arrays is 
partitioned into blocks corresponding to contributions from individual nodes, 
strain, or stress-resultant approximation functions. The expressions of the 
typical partitions (or blocks) are given in table Dl. Note that the order of the 


strain components in these partitions is e 


' 0 ' 


2e 


s6 9 


"0 9 


2k 


s0 ' 


2e 


s3 9 


and 2£0 3 . The order of the stress resultants is N s , Nq , n s q/ M g/ Mq , m s q/ 
and Qq. The order of the nodal displacement parameters is u, v, w, <J> S , and 


In table Dl, and Ar are the shape functions for both the stress re- 
sultants and strain components; are the shape functions for the generalized 

displacements; r is the radial coordinate; m is the number of displacement 
nodes in the element; [i] is the identity matrix; h is the number of param- 
eters used in approximating each of the stress resultants and strain components; 
(e) 

S is the element domain? 3 S = 3/3s; and 00 = 3/30. The range of the indices 

JL and j is 1 to 4 , and the range of the indices i and j is l_to m. The 
dots in the matrices refer to zero terms. The quantities and are 

defined in terms of the displacement components as follows: 


e, = 3 w 

1 R s 
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Table Dl.- Explicit Form of Typical Partitions of the Arrays [K], [K]> [r], [s], {q} 

and {I(H,X)} 























Table D1 . - Concluded 
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CO 
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CD 
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APPENDIX E - EVALUATION OF GLOBAL APPROXIMATION VECTORS 
BASED ON TWO-DIMENSIONAL SHELL FINITE E LE MENTS 


The global approximation vectors in equations (26) to (28) are obtained by suc- 
cessive differentiation of the governing finite-element equations (24) , with respect 
to the two parameters q and X, and solution of the resulting system of linear 
simultaneous algebraic equations. For individual finite elements, the recursion 
formulas for the global approximation vectors can be written in the following 
compact form: 


K 


9G, 


S + 


ax. 


Symmetric 


am ] 

ax! 


s m+n 


aq m dX n 


r E^ 


fp ( m+n ) ' 

H 

>-< 

^(m+n) 

< xj 




(El) 


where the range of I is 1 to 86 (6 is the number of parameters used in approxi- 

mating each of the stress resultants or strain components) and the range of I and 
J is 1 to 5m (m is the number of displacement nodes) . The total number of 
(m + n) combinations is K - 1, where Ji is the number of approximation vectors. 

The explicit forms of the first few values of are 

given in table El. 


In table El, a dot (*) over a symbol refers to a derivative with respect to q 
and a prime (') over a symbol refers to a derivative with respect to A. Note that 

the coefficient matrix on the left-hand side of equations (El) , which must be 

factored, is the same for each of the global approximation vectors. Hence, this 
matrix is factored only once regardless of the number of global approximation vectors 
generated. If the global approximation vectors are generated at A = q = 0, then 

{e} == 0, {h} = 0, and {x} = 0. Also, all the derivatives — ~{e}, — ;t{h}, and 

2 Yi 

— -{x} vanish, and equations (El) are thus simplified. The computational effort 
9A 

in evaluating the global approximation vectors is thereby reduced. 
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Table El.- First Explicit Forms of {P^ m+H ^}, {£) } , and (jj( m+n )} 
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SYMBOLS 


a 


radial distance for toroidal shell (see fig. 6) 


b o' b l' b 2 


parameters defining the elliptic profile of the tire cross section 
(see fig. 6) 


c. . ,d. . ,f . . 
i J il il 

c 44' C 45 ,C 55 


E ,E 
L' T 


shell stiffness coefficients (i,j = 1,2,6) 

transverse-shear stiffness coefficients of the shell 

elastic moduli of the individual layers in the direction of fibers 
and normal to it 


{e} 


vector of strain parameters for a shell element (see eqs. (11) 
and (24)) 


{£},{£} 


vectors of strain parameters for a one-dimensional shell element 
(see eqs. (13) and (14)) 


G rm ,G m m shear moduli in plane of fibers and normal to it 

LT TT * 

{©(X)}, vectors of nonlinear terms for a two-dimensional shell element 

(I!t(H,X)} (see eqs. (24)) 


{H} vector of stress-resultant parameters (see eqs. (11) and (24)) 

vectors of stress-resultant parameters for a one-dimensional shell 
element (see eqs. (13). and (14)) 


[k],[K],[K], 

[K],[K],[fC] 

W,[fe] 


total thickness of the shell 

elemental matrices (see eqs. (11) , (13) , (14), and (24)) 

linear stiffness matrices of the reduced system (see eqs. (21), 
(22) , (30) , and (31)) 


M s' M e' M se 

m 


m s ' m e 




W N se 

N 

W 


n 


bending (and twisting) stress resultants (see fig. 2) 
number of displacement nodes in the element 

intensity of external moments in the meridional and circumferential 
directions 

vector of nonlinear terms in the reduced system (see eqs. (32)) 
extensional stress resultants (see fig. 2) 

shape functions used for approximating the generalized displacements 

shape functions used for approximating the stress resultants and 
strain components 

Fourier harmonic in the circumferential direction 
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intensity of normal pressure loading 


p s ' p e' p 


intensity of external loading in the coordinate directions 
(see fig. 2) 


(q 2 ) 


vectors of normalized external loadings for individual elements 
(see eqs. (24) and (34)) 


q^,Qq transverse-shear stress resultants (see fig. 2) 

vectors of consistent external loadings and moment coefficients 
for a one-dimensional shell element (see eqs. (13) and (14)) 


q,q r q 2 

W 

R 1 ,R 2 

[r],[R],[R] 


loading parameters 

loading vector of the reduced system (see eqs. (23) and (29)) 

principal radii of curvature in the meridional and circumferential 
directions 

linear matrices of the shell element (see eqs. (11), (13), (14), 
and (24)) 


r 


normal distance from the shell axis to the reference surface 


h, number of global approximation vectors 

[s],[j?],[S] linear strain-displacement matrices of the shell element (see 
eqs. (11), (13), (14), and (24)) 


s 


meridional coordinate of the shell (see fig. 2) 


4 


number of parameters used in approximating each of the stress 
resultants and strain components within the individual elements 


U total strain energy of the shell 

u, v, w displacement components of the middle surface of the shell in the 

meridional, circumferential, and normal directions 


normal displacement at the center 


(X> 


vector of nodal displacements 


W,{X} 


vectors of generalized nodal displacement coefficients for a one- 
dimensional shell element 


x,y,z Cartesian coordinate system (see fig. 3) 

x^ coordinate normal to the shell middle surface 
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matrices of global approximation vectors defined in equations 
(16), (17), (18), (19), (26), (27), and (28) 


[r E ],[r H ],[r x ], 

[r £ ],[r H ],[r x ], 

[r E ],[r H ],[r x ] 


e s' e 0' 2e s0 

2£ s3 ,2e 03 


W 2lC s0 


LT 

5 

«)> 

V*0 

3 


e 


extensional strains of the middle surface of the shell 
transverse-shear strains of the shell 

circumferential (hoop) coordinate of the shell (see fig. 2) 
bending strains of the shell 

tracing parameter identifying all the anisotropic terms in the 
governing finite-element equations (see eqs. (13) , (14), and (24)) 

major Poisson's ratio of the individual layers 

dimensionless coordinate along the meridian (see fig. 3) 

rotation around the normal to the shell middle surface 

rotation components of the middle surface of the shell (see fig. 2) 

vector of amplitudes of global approximation vectors 

= 3/3s 

= 3/39 


Ranges of indices: 



1 

to 

a* 

I,J 

1 

to 

5m 

l,j 

1 

to 

6 


1 

to 

m 


Superscript t denotes transposition; subscript n 
harmonic. 


refers to the nth Fourier 
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Table 1.- Symmetry Conditions for Tires Exhibiting Rotational and/or 
Reflection Symmetry in Addition to Axial Symmetry 

[See fig. 3] 


Additional symmetries 
of tire and loading 

Tire type 

Symmetry conditions 

Rotational symmetry about 
an axis normal to the axis 
of revolution (e.g., 
z-axis) 

Ortho tropic 
or 

anisotropic 

At 5=0: 
u = v=<f> s = 4> e = o 
8 S - B e - 0 
2e s3 ' 2C 03 * ° 

Reflection (or mirror) 
symmetry with respect 
to x-z plane 

Ortho tropic 

At 5=0: 
u = <t> s = 0 

H s e = M se - 2 S - o 

2e s6 " 2K s0 = 2E s3 * ° 

Rotational symmetry 
about axis normal to 
axis of revolution and 
reflection symmetry with 
respect to x-z plane 

Ortho tropic 

i 

I 

i 

At -1 < 5 < 1: 

V = <t> 0 = 0 

N se " M se = Q e " 0 
2£ s0 = 2K s0 = 2£ 03 = ° 
At 5 = 0: 
u = <(> s = 0 
Q s = 0 

2£ s3 = ° 
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Table 2.- Symmetric and Antisymmetric Shell Variables 
[symmetric or antisymmetric with respect to 0 = o] 


Vector 

Symmetric variables 

Vector 

Antisymmetric variables 

{£} 

n 

e , , -2 e n , k , 

s,n 9 , n s0,n s,n 

K 0,n' 2K s0,n' 2G s3,n f 

- 2i 63,„ 

{E} 

n 

£ , , 2e , k , 

s,n 0,n s0,n s,n 

^0 ,n' 2lC s0 ,n' 2 ^s3,n' 

2e 03,n 

{H} 

n 

N s,n' N 0,n' ^s 0 , n ' M s,n' 
M 0,n' ^s0,n' ^s,n' ^0,n 

i 

N c N , N , M , 

s , n b So r n s ; n 

Wq F M A F Q F Qa 

0 ,n s0 ,n s ,n 0 / n 

{X} 

n 

U , V , W , <J> , $ 

n n n s / n 0 ^ n 

{X} 

n 

U n' V n' W n' ^s,n' ^0,n 

<«> n 

s f n v ^ n n s / n b ^ n 

{ * } n 

P s,n' P 0,n' ^n' ™s,n' m 0,n 


Table 3.- Symmetric and Antisymmetric Tire Quantities 

[n = 0] 


Vector 

Symmetric variables 

Vector 

Antisymmetric variables 

«>o 

e s,o' e e,o' K s,o' K e,o' 
2£ s3,0 

fE} 0 

2e s0,O' 2lC s0,O' 2G 03,O 

{H} 0 

”s,0' N 9,0- M s, 0' M e,o' 
®s,0 

™»o 

1 

N s0,O' M s0,O' Q 0,O 

{X} 0 

"o' “o' *5,0 

{x) o 

"V -*0,0 

<«>0 

P s,0' P 0' m s,0 

^>o 

- p 0,o' “ m 0,o 
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Table 5.- Accuracy and Convergence of Solutions Obtained With Proposed Strategy 
for Orthotropic Tire Subjected to Uniform Inflation Pressure and 

Unsymmetric Localized Loading 


See fig. 12; for uniform inflation pressure only (q 1 = q 2 = 0) , 
w w, 

= 0.08935 and — ~ = 0.3633 
h h „ , 3 


Number of 

approximation vectors 

* 

w 

XjJ * 

10x^ 

h 

U 

c 

h~ 

v 3 

3 

-0.6430 

0.8820 

0.3464 

6 

-.6499 

.8907 

.3572 

10 

-.6548 

.8896 

.3570 

15 

-.6569 

.8842 

.3568 

21 

-.6554 

.8917 

.3569 

Direct finite-element 
solution 

-0.6553 

0.8918 

0.3569 


* 

See figure 3. 



Tread 



Figure 1.- Components of an aircraft tire. 











Figure 3.- Axes of rotational symmetry and plane-of-ref lection 

symmetry for tires. 




Analysis model 


Typical meridian 


Figure 4.- Axial symmetry for orthotropic and anisotropic tires. 
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(a) Symmetric generalized 
displacements . 


(b) Antisymmetric generalized 
displacements . 


Figure 9.- Accuracy of solutions obtained with proposed strategy for two-layered 
anisotropic tire subjected to normal loading and modeled by semianalytic 
finite elements and n = 2. p = p Q cos 20. 
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Figure 11.- Convergence of maximum normal displacement w c and total 
strain energy U obtained with proposed strategy and with Taylor 
series expansion for two- layered anisotropic tire subjected to 
localized loading. 







Figure 12.- Concluded 
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